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ASCA Discovery of a Be X-Ray Pulsar in the SMC: 

AX J0051-733 
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Abstract 

ASCA observed the central region of the Small Magellanic Cloud, and found a hard X-ray source, AX 
J0051— 733, at the position of the ROSAT source RX J0050.8— 7316, which has an optical counterpart of a 
Be star. Coherent X-ray pulsations of 323.1 ± 0.3 s were discovered from AX J0051— 733. The pulse profile 
shows several sub-peaks in the soft (0.7-2.0 keV) X-ray band, but becomes nearly sinusoidal in the harder 
(2.0-7.0 keV) X-ray band. The X-ray spectrum was found to be hard, and is well fitted by a power-law 
model with a photon index of 1.0 ± 0.4. The long-term flux history was examined with the archival data 
of Einstein observatory and ROSAT; a flux variability with a factor > 10 was found. 
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1. Introduction 

Systematic X-ray source surveys in the Small Magel- 
lanic Cloud (SMC) started with the Einstein Imaging 
Proportional Counter (IPC: see e.g. Wang, Wu 1992), 
followed by the ROSAT Position Sensitive Proportional 
Counter (PSPC: Kahabka, Pietsch 1996; Kahabka et al. 
1999). Optical identifications, however, require more ac- 
curate determinations of the X-ray source positions, and 
have been made by Cowley et al. (1997) and Schmidtkc 
et al. (1999) with the ROSAT High-Resolution Imager 
(HRI). They conducted optical photometry and spec- 
troscopy studies with the Cerro Tololo Interamerican Ob- 
servatory (CTIO) and identified X-ray sources with opti- 
cal counterparts. Most of these were proven to be super- 
nova remnants (SNRs) in the SMC, or foreground bright 
stars in our Galaxy, while four were found to have Be 
star counterparts within their error regions. Hughes and 
Smith (1994) also made an optical counterpart search for 
selected X-ray sources, and found two Be star candidates 
at the X-ray positions. 

X-ray binary pulsars with a Be star companion (here 
Be-XBPs) are thought to comprise the majority of X-ray 
binary pulsars (XBPs). However, possibly due to their 
transient nature, most of the Be-XBPs have not yet been 
discovered in the past X-ray observations. The SMC re- 
gion would not be an exception. Since a pulsation search 
from X-ray sources with a Be star companion is a direct 
approach to discover Be-XBPs, we have performed tim- 
ing studies for X-ray sources with a Be star companion 
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reported by Cowley et al. (1997), Schmidtke et al. (1999) 
and Hughes, Smith (1994). The studies were made using 
archival data of the Japanese X-ray satellite ASCA. 

In this Letter, we report on the discovery of X-ray pul- 
sations from AX J0051-733 = RX J0050.8-7316, which 
was the brightest X-ray source among the optically iden- 
tified Be stars. We also report on the long-term flux vari- 
ability from Einstein, ROSAT and ASCA observations. 

2. Observation 

AX J0051-733 was included in the ASCA obser- 
vation of the radio supernova remnant N19, on 1997 
November 13-14, pointing at the direction of (00 h 47 m 16 s , 
-73°08'30") (here and after, we use J2000 coordinates). 

ASCA (Tanaka et al. 1994) is equipped with four iden- 
tical X-ray telescopes (XRTs: Serlemitsos et al. 1995) 
sensitive to photons with energies 0.4-10 keV. Four in- 
dependent detectors, the two Solid-state Imaging Spec- 
trometers (SISs: Burke et al. 1991) and the two Gas 
Imaging Spectrometers (GISs: Ohashi et al. 1996; Mak- 
ishima et al. 1996), are placed at the foci of XRTs and 
are operated in parallel; hence, four independent data 
sets can be provided. However, we analyzed only the 
GIS data, because AX J005 1 — 733 was outside of the SIS 
field of view. 

Each of the GIS detectors was operated in the PH 
mode with the standard bit-assignment that provides a 
time resolution of 62.5 ms and 0.5 s in high and medium 
bit-rate, respectively. Data taken at a geomagnetic cutoff 
rigidity lower than 4 GV, at an elevation angle less than 
5° from the Earth and during passage through the South 
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Atlantic Anomaly were rejected. The particle events were 
also removed by a rise-time discrimination method. Af- 
ter these filterings, the total available exposure time of 
each GIS was w 43 ks. 
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3. Analysis and Results 

3.1. X-Ray Images and Source Identification 

We made GIS2 and 3 images of the N19 region, and 
found a bright point source at an off-axis angle of ~ 16' 
in the GIS fields. The position was determined to be 
(00 h 50 m 50 s , -73°16'04"), with a possible error of - T. 
The larger position error than the nominal value of ~ 40" 
(Makishima et al. 1996) is due to its large off-axis angle. 
Within this error region of AX J0051— 733, we found a 
ROSAT source, RX J0050.8-7316, with coordinates of 
(00 h 50 m 45?3, -73°15'54"), which was optically identified 
to be a Be star of V ~ 15.4 mag, and B - V = -0.04 
mag (Cowley et al. 1997). 



3.2. Timing Analysis 

Since the shape of the point-spread function is elon- 
gated at an off-axis position, like AX J0051— 733, we ex- 
tracted photon events for GIS 2 and 3, each from an 
elliptical region with major and minor axes of ~ 5' and 
~ 2', respectively. 

After a barycentric photon arrival time correction, we 
performed FFT (Fast Fourier Transformation) analyses; 
we found a peak in the power spectrum. To maximize the 
signal-to-noise ratio of this peak, we tried FFT analyses 
with many trial energy bands. Figure 1 shows the best 
power spectrum using the 1.5-6.0 keV energy band. A 
total of 541 (GIS 2) and 488 (GIS 3) photons were used 
for this FFT analysis. 

A clear peak is noted at a frequency of 3.095xl0~ 3 
Hz, corresponding to a pulsation period of ~ 323 s. The 
chance probability to detect such a strong power signal at 
any frequency given in figure 1 is only ~ 10~ 10 ; hence, we 
can safely conclude that the peak is not due to Poisson 
fluctuations, but is a real signature. 

To estimate the pulsation period more precisely, the 
epoch-folding technique was applied, assuming that no 
period change occurred during the observation. The most 
likely period was determined to be P — 323.1 ± 0.3 s. 
Figure 2 shows the folded light curves in the soft (0.7-2.0 
keV: upper panel) and hard (2.0-7.0 keV: lower panel) X- 
ray bands. The pulse profiles are different between these 
two bands; we can see three or four sub-peaks in the 
soft band, while a quasi-sinusoidal curve can bee seen in 
the hard band. The pulse fractions, defined as (Pulsed 
intensity) /(Total intensity), are ~ 56% and ~ 80% in the 
soft and hard bands, respectively. 
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Fig. 2.. Folded light curves of AX J0051-733 in the 0.7-2.0 
keV (upper panel) and 2.0—7.0 keV (lower panel) bands. The 
background levels are represented by the dotted lines. 
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Fig. 3.. Upper panel: Phasc-avcragcd energy spectrum of 
AX J0051-733 obtained from GIS2. The solid line shows the 
best-fit power-law model. Lower panel: Each residual from 
the best-fit model. 



3.3. Spectral Analysis 

Unfortunately, AX J005 1—733 is located only ~ 6' 
away from the calibration isotope of GIS 3. We therefore 
did not use the GIS 3 data for the spectral study in order 
to avoid possible contamination from the isotope. We ex- 
tracted a source spectrum of GIS 2 from the same region 
as that used for the timing analysis and selected a section 
of nearby sky as a background region. Figure 3 shows 
a phase-averaged source spectrum. Since the spectrum 
showed neither an emission line nor any other prominent 
structure, we applied a simple power-law model with in- 
terstellar absorption for the fitting. The best-fit model 
and parameters are given in figure 3 and table 1, respec- 
tively. 

To determine any possible spectral changes during the 
pulse, we divided the data into two phases, phase 0.0- 
0.5 (high state) and 0.5-1.0 (low state), and fitted two 
separate spectra. Since the statistics of these spectra 
were limited, we fitted them to a power-law model fixing 
the absorption column density to the best-fit value of the 
phase-averaged spectrum. The best-fit parameters are 
summarized in table 1. The large errors make it unclear 
whether or not the photon index changed with the pulse 
phase. 

4. Long-Term Flux History 

To investigate any long-term flux variability, we ac- 
cessed the HEASARC archive system for the Einstein 
and ROSAT observations, and found that AX J0051-733 
was included in 1 Einstein IPC, 9 ROSAT HRI and 6 
ROSAT PSPC observations. For each observation, we 
extracted the source and background events from a cir- 
cular region of 1' radius and an outer region of I'—l'A 
radius, respectively. After background subtraction, we 



Fig. 4.. Flux variation history of AX J0051-733. Data are 
obtained by the Einstein IPC (square), ROSAT PSPC (cir- 
cle), ROSAT HRI (triangle), and ASCA GIS (star). The data 
point given by the arrow (MJD = 49082.2-49104.0) is the 
ROSAT PSPC observation that showed no excess X-ray from 
AX J0051— 733, hence we can only determine the upper limit 
of the flux. 



estimated the source flux using PIMMS software, assum- 
ing that the spectral shapes were the same as that ob- 
tained with ASCA. 

In figure 4, we summarize the X-ray fluxes from AX 
J0051-733. The ASCA flux at MJD = 50765.2-50766.3 
was higher than those of Einstein obtained at MJD = 
44189.1-44189.9 and those of ROSAT at 49082.2-49104.0 
by a factor > 10. Therefore, we conclude that AX 
J0051 — 733 has been highly variable, although a possibil- 
ity that low flux periods are due to eclipses is not ruled 
out. 

5. Discussion 

5.1. Is AX J0051-733 a Be Star Pulsar ? 

We have discovered coherent X-ray pulsations from a 
Be star binary AX J0051-733 = RX J0050.8-7316. The 
long pulsation period (~ 323 s), large flux variation in the 
long time span, flat power-law spectrum (photon index ~ 
1) and existence of a Be star companion strongly support 
that AX J0051-733 is a Bc-XBP (Nagase 1989). 

In the long-term history, we found 7 X-ray outbursts 
with flux > 2 x 10~ 12 erg s _1 cm" 2 . Most of the known 
Be-XBPs exhibit highly eccentric orbits, and flaring ac- 
tivities are often taken place when a neutron star passes 
near the peri-astron or across the equatorial plane, where 
the circumstellar gas density is expected to be high. We 
suspect that the occasional flares from AX J005 1—733 
are related to a possible eccentric orbital period. How- 
ever, no clear periodicity is noted in these flare intervals. 
We, nevertheless, folded the long-term light curve with 
a period of 185 days, which is near the expected orbital 
period based on Corbet's empirical relation between the 
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orbital and pulse periods in a Bc-XBP (Corbet 1984). 
Although a hint of clustering of the outburst around an 
interval of about 185 d was found, as is given in figure 
5, no definite argument for an orbital period was ob- 
tained. To establish the Be-star pulsar binary scenario 
and its structure, an orbital period determination is es- 
sential. We encourage further monitoring observations of 
AX J0051-733. 

5.2. Pulsation Search from the Other Be Star Binary 
Candidates 

The discovery of coherent pulsations from AX 
J0051— 733, which was optically identified to be a Be 
star binary in the SMC, leads us to suspect that the 
other Be star binary candidates (RX J0051.9— 7311, RX 
J0052.9-7158, and RX J0058. 2-7231: sec Schmidtke 
et al. 1999; RX J0103.2-7209 and RX J0106. 2-7205: 
see Hughes, Smith 1994) also exhibit coherent pulsations 
(Be-XBPs). We thus extended our pulsation search to 
these sources. 

RX J0051.9-7311 and RX J0058. 2-7231 were found 
in the ASCA archive; both are found to have hard X-ray 
spectra (Yokogawa et al. 1999). Since the hard X-ray 
spectra are one of the signatures of XBPs, we extensively 
searched for pulsations, but found none, mainly due to 
their limited statistics. In fact, the total photon num- 
bers for RX J0051.9-7311 and RX J0058. 2-7231 in the 
1.5-6.0 keV band are respectively 361 and 273 (including 
background counts), which are only ~ 1/3 and ~ 1/4 of 
that of AX J0051-733. 

RX J0103. 2-7209 = 1SAX J0103. 2-7209 is located in 
the radio source SNR 0101—724 and has been found to 
exhibit coherent pulsations of 345.2 s with BcppoSAX 
(Israel et al. 1998). This source was seen in the three 
ASCA archival data sets. Coherent pulsations of 348.9 
s were found in the data set obtained on 1996 May 21, 
the highest statistics from this source among the three 
observations (Yokogawa, Koyama 1998). 

RX J0106.2— 7205 showed no detectable X-ray emis- 
sion in two ASCA observations, while RX J0052. 9-7158 
has not been observed with ASCA. 

In summary the brightest two X-ray sources exhibit 
X-ray pulsations, while three weaker sources show no de- 
tectable pulsations. This latter fact is most probably due 
to their limited statistics. The remaining source was not 
observed with ASCA. Since Be-XBPs often exhibit X-ray 
outbursts, we may have a chance to detect X-ray pulsa- 
tions even from the three weaker sources during their out- 
burst episodes. Thus, long-term monitoring over a wide 
area of the SMC would be very fruitful in establishing 
the statistics of the Be-XBP population in a neighboring 
galaxy. 
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Fig. 5.. Same as figure 4, but the data are folded with a trial 
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Table 1. Best-fit parameters to the spectra of AX J0051 — 733. 



Phase 


Photon index* 


Column density* 


Flux 1 " 


Luminosity'* 


X 2 /d-o.f. 






(H cm" 2 ) 


(erg s _1 cm~ 2 ) 


(erg s" 1 ) 




Average 


1 +u ' 4 


4.4+^ xlO- 


3.4xl0 -12 


1.6x10* 


36.04/33 


(0.0-1.0) 












High 


i-ii° J 


4.4 xlO 21 (fix) 


4.3xl0~ 12 


2.0 xlO 36 


31.37/19 


(0.0-0.5) 












Low 


0.8i ( ^ 


4.4 xlO 21 (fix) 


2.5xl0~ 12 


1.2 xlO 36 


12.52/12 


(0.5-1.0) 













*: Errors are at 90 % confidence level. 

f: In the 0.7-10.0 keV band. 

J: Absorption-corrected value at the source distance of 60 kpc (Mathewson et al. 1988). 



